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Abstract 

We present a study of connectivity percolation in suspensions of hard spherocylinders by means of 
Monte Carlo simulation and connectedness percolation theory. We focus attention on polydispersity in 
the length, the diameter and the connectedness criterion, and invoke bimodal, Gaussian and Weibull 
distributions for these. The main finding from our simulations is that the percolation threshold shows 
quasi universal behaviour, i.e., to a good approximation it depends only on certain cumulants of the 
full size and connectivity distribution. Our connectedness percolation theory hinges on a Lee-Parsons 
type of closure recently put forward that improves upon the often-used second virial approximation [T|. 
The theory predicts exact universality. Theory and simulation agree quantitatively for aspect ratios in 
excess of 20, if we include the connectivity range in our definition of the aspect ratio of the particles. 
We further discuss the mechanism of cluster growth that, remarkably, differs between systems that are 
polydisperse in length and in width, and exhibits non-universal aspects. 
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1 Introduction 


Composite nanomaterials have long attracted attention because of their potential application for in¬ 
stance in electronics, display technology and photovoltaics [2]. Of particular interest in this context are 
their heat and charge transport properties [3]. Adding a sufficient amount of electrically conductive 
nanoparticles, such as carbon nanotubes or graphene, to an insulating polymer matrix produces a 
conductive composite the conductivity of which can be tuned by the choice of filler type, filler loading 
and processing pf|. For many technological applications, the minimum filler loading required to reach 
a conductive state, the so-called percolation threshold, is desired to be as low as possible [5]. Rod-like 
particles are particularly suitable for this kind of application since they present very low percolation 
thresholds (6] • For the purpose of the rational design of such materials it is crucial to be able to describe 
and predict the percolation threshold of assemblies of filler particles, and understand the underlying 
mechanisms of the buildup of the system-spanning network required for effective conduction. 

In experimental reality the properties of the filler nanoparticles are not always well controlled [7j. 
Indeed, they are usually chemically and otherwise polydisperse, that is, consist of a mixture of particles 
of different dimensions and conductive properties. This complexity makes prediction of the percolation 
threshold and of the network structure very difficult, not least because of the huge parameter space. In 
this paper we present a simulation and theoretical study of percolation in dispersions of polydisperse 
nanorods, specifically allowing for hard core interactions and targeting aspect ratios that are of an 
intermediate range, i.e., not in the scaling limit mmm- 

Even though we find qualitative agreement with work on polydisperse ideal (penetrable) rods [10| 
and very long hard rods [81 [9] (showing that the percolation threshold obeys laws that within a good 
approximation depend only on a few moments of the full distributions functions) quantitatively our 
results are very different. In fact, we find strong deviations in the dependence of the percolation 
threshold on the appropriate measures for the mean aspect ratio and connectivity of the particles. 
Finally, we find that the network connectivity properties are affected differently by variabilities in 
length, diameter and connectivity criterion. 

It is important to point out that the model systems that have been studied in the literature so far 
usually capture only one or a few aspects relevant to experimental reality. A very large focus is on the 
particle shape, attractive interactions and aspect ratio m ■ While there is a huge body of literature 
dealing with monodisperse systems, relatively little attention has been paid to polydisperse systems 
[m uni nn usi uni im nn nm eoi eu [221 o eh 123 hu eh ez]- Recently, Chatterjee [28] and Otten 
and Van der Schoot M have developed theories of continuum percolation that take polydispersity 
into account and predict universal scaling laws for the percolation threshold. These predictions have 
only to a small extent been tested numerically. 

In a recent simulation study, Nigro and co-workers confirm that for hard and penetrable rods that 
are polydisperse only in length the percolation threshold depends only weakly on the exact shape of 
the length distribution m ■ a similar finding was obtained by Mutiso and collaborators for mutually 
penetrable length and width polydisperse rods [29]. They also find that finite-aspect-ratio corrections 
on the predictions of Otten and Van der Schoot are quite significant up to aspect ratios of about one 
hundred mm- 

Here we go considerably beyond the scope of earlier work, and report on simulation results for three 
different types of polydispersity that we investigate separately. The coupling between different kinds of 
polydispersity, predicted to be relevant for many experimental systems |9], is postponed to future work. 
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Figure 1: Definition of the particle dimensions and connectivity range. 

We show that the different kinds of polydispersity exhibit non-trivial universal behaviour. We invoke 
a treatment of connectedness percolation theory of hard rods recently put forward by us [Tj , which 
is aimed at predicting finite-aspect-ratio corrections rather than obtaining them phenomenologically 
front simulations, as was done in [29j . 

The remainder of this paper is arranged as follows. We present in section 2 the methods imple¬ 
mented in the Monte-Carlo simulations. Section 3 deals with the derivation of our version of connec¬ 
tivity percolation theory for polydisperse spherocylinders, including the Lee-Parsons approximation. 
Finally, we focus in section 4 on both numerical and theoretical results, first about the percolation 
thresholds and then about the cluster mechanisms. We end the paper with conclusions and a summary 
of the main findings in section 5. 

2 Simulation Methods 

We consider hard spherocylinders consisting of cylinders of length Li and diameter Di , each capped by 
two hemispheres of the same diameter. See fig. [0 The particles are not allowed to overlap, but do not 
directly interact with each other when they are not in contact. The corresponding interaction potential 
is therefore either zero or infinite, making their resulting equilibrium properties temperature indepen¬ 
dent. We initialize a simulation box in which around 10, 000 spherocylinders are perfectly aligned and 
regularly placed on square lattices spaced from each other along the rod direction. At each simulation 
step, the particles are then randomly rotated and translated. Equilibration is monitored by computing 
the nematic order parameter, which is expected to reach a constant value at the equilibrium (0 in the 
isotropic phase). Once the system is equilibrated, we generate ca. 5000 independent configurations of 
the system and average all quantities of interest over those configurations. In order to detect overlap¬ 
ping particles efficiently, the box is divided into a fine grid ED], where the unit cell length is chosen 
equal to the greatest rod diameter in the system, so that the computational cost increases linearly 
with the number of particles. This method is very fast but rather expensive in terms of memory. 

We take polydispersity into account by assigning to each rod a length and a diameter according 
to a probability density function Vl and Vd- To generate a finite number of rods from a continuous 
distribution, we define an interval Ax where x = L, D stands for length and width, and force the 
system to contain NV x (x o)Ax rods, whose dimension lies between xq and xq + Ax. This method 
turns out to give much more accurate results than simply drawing the rod dimensions directly from 


3 








Figure 2: Snapshot of an equilibrated configuration of diameter polydisperse spherocylinders at the critical volume 
fraction (left) and of the largest cluster within this particular configuration (right). Distribution is of the Weibull form 

P(D) oc with a = 2.83 and /3 = 1.16. Lengths and connectendess distances are all fixed to L = 15 

and £ = 0.2. 


the distribution under study. In this work, the mean aspect ratio lies around L/D = 15 and the widest 
distributions we considered spread up to an aspect ratio of approximatively L/D = 80 for the very 
longest rods. In order to clearly distinguish between the effects of length and diameter polydispersity 
we choose only uncorrelated distributions. As already advertised, this assumption does not necessarily 
apply to all experimental systems. 

Connectedness percolation requires the definition of an inter-particle connectedness criterion. We 
define for each rod i a spherocylindrical shell of length Li and diameter Di + that contains the 
particle, where the connectivity parameter obeys some distribution function V(£/). Two particles 
are then connected if their surrounding shells overlap. Clusters are defined by contiguous pairwise 
connections. We define a configuration percolating if one of its clusters is connected to its image 
under periodic boundary conditions. To every configuration corresponds a percolation probability 
that is either 1 (it percolates) or 0 (it does not percolate), and averaging over many configurations 
we compute a global continuous percolation probability for a particular system. A typical snapshot of 
such a sample and of its corresponding largest cluster is shown in fig. [2] 

In order to estimate the percolation threshold for a particular length, diameter or connectivity 
distribution, we perform simulations using this distribution for a range of rod volume fractions. The 
volume fraction </> is defined with respect to the hard core volume of the particles and does not take 

N 

into account the connectivity shell: <fi = Y/ v i where Vi = jLiDf + ^D\ is the volume of the 
particle i and a is the simulation box length. The percolation probability in a finite system is a 
sigmoidal function of the volume fraction running from 0 to 1. Its transition steepness increases with 
the box size and reaches a Heaviside step function in the limit of infinite box volume. The curves that 
correspond to different box volumes cross each other slightly below the concentration at which the 
probability reaches the value of 0.5. As we are interested in the scaling behaviour of the percolation 
threshold with the aspect ratio and cumulants of the size distribution of the particles, we do not need 
very accurate estimates. Hence, we ignore finite size effects and assume that the percolation threshold 
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is the volume fraction corresponding to a percolation probability of 0.5. We verify that our box is 
sufficiently large to ensure that the percolation probability goes from 0.2 to 0.8 within a maximal 
volume fraction range A 4> max = 0.005. We assume this criterion to be restrictive enough for our 
results to achieve a satisfactory accurracy. 

3 Theory 

Percolation of clusters of nanoparticles in a fluid background medium can be investigated theoretically 
invoking what in essence is liquid state integral equation theory [5j. The theoretical framework is called 
connectedness percolation theory and it has been applied to hard and soft rod-like particles mmm- 
Here, we follow the same recipe, except that we will not rely on the second virial approximation that 
becomes exact in the limit of infinite aspect ratio. Instead we opt for a closure that is was recently 
shown to provide an accurate description of percolation of monodisperse, hard rods with an aspect 
ratio larger than roughly ten [1]. 

Within the framework of connectedness percolation theory, the cluster size S can be expressed 
in terms of a function T, averaged over all of the attributes of orientation vector u and dimensions 
x = (L, D , £) of the particles [9], 

5=(T(x,u)) x>u . (1) 

The function T itself is a pair connectedness function averaged over its attributes, and the solution of 
a generalised connectedness Ornstein-Zernike equation 

T(x,u) - p/c' + (0,x,x',u')T(x',u')) =1 • (2) 

Here, p is the number density of particles and C + = (7 + (q, u, u') the spatial Fourier transform of the 
connectedness direct correlation function, that is, the direct correlation function for particles that are 
part of the same cluster, and q the wave vector. To average over the entire volume of the system, we 
have to take the zero wave vector limit, q —> 0. 

In the isotropic phase the rods are randomly oriented, implying that T(x, u) = T(x), which in 
turn allows us to redefine C + as its average over the possible orientations, producing the simplified 
connectedness Ornstein-Zernike equation 

T(x)-p(c + (0,x,x')T(x')) x/ =1 (3) 

This equation needs to be closed and we follow Schilling et al. [1J by invoking the following Ansatz: 
C' + (0,x,x') = r(»/+(0,x,x') pp. Here, f + = e~^ u+ is the connectedness Mayer function and /+ its 
spatial Fourier tranform |9j , with [3 the reciprocal thermal energy and u + the so-called connectedness 
potential, and r(</>) a coefficient that depends on the volume fraction </> of the particles. The functional 
form of r is obtained from the Lee-Parsons expression for the excess free energy, interpolating between 
the Percus-Yevick equation of state for hard spheres and the Onsager equation of state for hard rods 

[3U [3211m EH [35]. Within this Ansatz , we have 

r « - (iv|p < 4 > 
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Notice that in the limit </> —> 0, T —> 1 and we obtain the second virial theory that is valid in the 
Onsager limit of very slender rods. As we demonstrated recently [1], corrections to the Onsager limit 
are significant for aspect ratios below a few hundred. 

For the case of hard rods, the connectedness Mayer function is 1 if the distance between two 
rods is between D and A = D + £ and 0 otherwise. The zero-wave vector Fourier transform of 
the connectedness Mayer function / + (0,x, x') can be separated into contributions from interactions 
between the different portions of the spherocylindrical particles. We use Onsager’s expression for the 
excluded volume of two hard spherocylinders to obtain [3B] 


/+(0, x, x') = LL'^-^-fa + (L + L') 


A +A' 




fw 


fio — 


A +A' 


D + D' 


foo 


foo 


(5) 


where A = D + £ and the coefficients /oo = 4-7r/3, /io = 7r and /n = n/2 respectively indicate the 
cylinder-cylinder, cylinder-hemisphere and hemisphere-hemisphere contributions. Note that because 
u + is infinite for particles with overlapping hard cores and for those that are not connected, and zero 
for connected ones, /+ is essentially the difference between the excluded volume of two particles with 
hard-core radius A minus that of particles with hard-core radius D. 

Let us first focus on length polydispersity alone, and set D = D' and A = A'. In this particular 
case eq. [5] becomes 

f + (L, L ') = LL'ivifu + (L + Z/)w 2 /io + W 3 / 00 , (6) 

where = A” — D n are differences between powers of the interaction ranges A and D. Inserting this 
into equation [3j we find 


T{L) - pT^XL^/n + oj 2 f 10 ) (LT(L)) l - pT^L^f 10 + W 3 / 00 ) {T(L)) l = 1. (7) 


The two unknown coupled quantities are (T(L)) l and (LT(L)) L , therefore we need two independent 
equations relating them. The first one we obtain by averaging equation [71 whereas the second we 
derive by first multiplying equation [7] with L and averaging the resulting equation. As these relations 
are linear, we can summarize them by defining two vectors X_ and Y_ and a matrix M : 


X = 

(T(L)) 

, Z = 

1 

, M = 

1 — kai —kj3i 


_ ( lt(l )} . 




—ka 2 1 — k/3 2 


MX = Y (8) 


where we use the notation k = pT(</>), a\ = (L)u> 2 fio + ^ 3 / 00 , a 2 = ( L 2 )w 2 /io + ( L)u} 3 f 00 , /3i = 
(L)wi/n +W 2 / 10 , and /3 2 = (L 2 )wi/n + (L)oj 2 fio- The cluster size S = ( T(L )} is the first element 
of X — M _1 F. Each element of the matrix M _1 is a fraction whose denominator is the determinant 
of M- 

Therefore, S can be written as 


S = 


E A npq (L n )DP £* 

n ,p,q 


det M 


(9) 


where A npq are coefficients which involve the quantities k and fij. However, their exact expressions 
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are not important for the remaining calculations, so we do not reproduce them here. The percolation 
threshold is the volume fraction <j) c for which S diverges, or, equivalently, that makes detM = 1 — 
k(a\ + P2) + k 2 (aif82 — a2P1) vanish. Therefore, we determine <f> c by solving this simple second-order 
polynomial equation for k and by then relating k and <j>, <pT((f>) = k (v), where (v) = j ((L) D 2 + | D 3 ) 
is the average volume of a particle. The polynomial equation yields two solutions k±. We define then 
7± = &± ( v )• 

The final equation that we have to solve is therefore 

^7± + ^ 4 % ~ (! + 27±)<^>c + 7 ± = 0 , ( 10 ) 


which has two solutions for (f> c , 




1 + 27 ± ± t/T + 7± 

27± + | 


( 11 ) 


We have four solutions for <j> c but we only keep the positive solution below unity for obvious physical 
reasons. 

The analysis of the cases of diameter and connectedness polydispersity are completely analogous 
to that of length polydispersity, that is, the same method is applied starting from eq. [5] setting L = L' 
and ? = ?', and L = L' and D = D', respectively. M then turns into a 3 x 3 matrix for diameter 
polydispersity and a 4 x 4 one for connectedness polydispersity. This leads to third and fourth order 
polynomial equations that need to be solved. It is important to note that the percolation threshold 
will then explicitly depend on the higher order moments ( D n ) with n < 4 and (? m ) with m < 6, 
respectively. 


4 Results 


4.1 Percolation thresholds 

If we invoke the second-virial approximation and neglect the end-cylinder and end-end interactions, 
the percolation threshold, 4> C: becomes proportional to the reciprocal weight average length (T)” 1 = 
( L ) /^L 2 ), the mean square width (H 2 ) and a measure for the mean reciprocal connectivity length 

(?) + V(? 2 ) > depending on the type of polydispersity mm 0;. Although this approximation 
produces results that are not very accurate for aspect ratios that are not huge, it is useful to take 
it as a reference because it shows what cumulants of the full distribution are expected to govern the 
percolation threshold. We note in this context that the nature of length and diameter polydispersity 
is fundamentally different from that of connectedness distance polydispersity. The first two relate to 
polydispersity in the particle dimensions and hence in the interactions between particles, whilst the 
third one is a polydispersity in the electrical connectivity length scale only, which does not affect the 
structure of the liquid but only the resulting cluster size distribution. 


To separate these various effects, we define two new dimensionless quantities, being \ — (L) w / \J{D 2 ) 
and A = (?) + \J (? 2 ) /2 yj ( D 2 ). The former, y, becomes equal to either (L) w /D or L/^Wj de¬ 
pending on the type of polydispersity and is analogous to the aspect ratio, as we can not define a 
unique aspect ratio in polydisperse systems. The latter, A, becomes (?) + (? 2 ) /2 D for rods with a 

monodisperse diameter, and represents a characteristic connectedness shell thickness compared to the 
particle diameter. 
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Figure 3: Percolation threshold <f> c for length polydispersity as a function of the reciprocal aspect ratio 1 = D/ (J■ 
Results are indicated from simulations (full dots), Lee-Parsons theory as well as the second-virial approximation (empty 
dots) for £/D = 0.2 and various distributions V(L). A remarkable universal scaling with x~ ! is observed in the three 
cases. See also the main text. 


We compare the percolation threshold obtained from our simulations, prediction based on the 
theory presented in the previous section and those from the second-virial theory for which T = 1, as 
a function of \ l f° r length and diameter polydispersity and as a function of A” 1 for connectedness 
polydispersity. 

For all three kinds of polydispersity, we tested bidisperse, Gaussian and Weibull distributions. The 
first describes binary mixtures, the second seems relevant as Gaussian distributions are common in 
many fields of physics, and the third has been experimentally observed in polymer-fiber composites 
that are polydisperse in length [87l 138] , 

Fig. [3] shows the percolation threshold for the case of length polydispersity as a function of the 
inverse aspect ratio x • Both theory and simulations display a remarkable universal behaviour of 
the percolation threshold as a function of (L) . The three different distributions that we tested are 
very different in shape, but we find that the percolation threshold is not sensitive to this if expressed 
in terms of the weight average length (L) w . This was also shown by Nigro et al. [TO] for penetrable 
particles and in a more limited fashion for hard particles. We find that the Lee-Parsons theory and our 
simulation results quantitatively converge in the range \ > 20 and are in qualitative agreement below 
that, whereas the theoretical prediction derived from the second virial approximation deviates notably 
even for relatively large aspect ratios. For monodisperse rods, the Lee-Parsons approach produces 
quantitative results already for y > 10 [T]. In polydisperse systems, the discrepancies between theory 
and simulations extend to larger average aspect ratios because of the shorter rods that are also present 
in the system. Note that the simulation results are always below the theoretical prediction. Hence, 
composite materials that contain fibres of short aspect ratio do not need as high filler loadings as 
theoretically expected in order to become conductive. 

Also for diameter polydispersity the second order cumulant is expected to be the most relevant one 
according to the second-virial theory, at least if we neglect end effects [9]. This is confirmed in fig. [4] by 
both our simulation results and the more accurate predictions of Lee-Parsons theory albeit that we do 
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Figure 4: Percolation threshold <fi c for diameter polydispersity as a function of x ~ 1 = y ( D‘ 2 )/L . Results are indicated 
from our simulations (full dots), Lee-Parsons theory as well as the second-virial approximation (empty dots) for L/£ = 75 
and various distributions V(D). The scaling with 1 is not anymore universal, higher order cumulants matter. See 
also the main text. 
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Figure 5: Percolation threshold (f) c for polydispersity in the connectedness distance £, as a function of A 1 = 
2D/ |"(£) + yj (^ 2 )j. Results from simulations (full dots), Lee-Parsons theory as well as second-virial approximation 
(empty dots) for L/D = 15 and various distributions 'P(£). A universal scaling with A -1 is also obtained. See also the 
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not observe strictly universal behaviour. This suggests that higher order cumulants must be important 
too, at least for aspect ratios below 25. Remarkably, even the second-virial prediction yields results 
that are in almost quantitative agreement with the simulations for \ > 15. The Lee-Parsons theory is 
in very good agreement with the simulation results even for relatively short particles. We note that 
this theory underestimates the percolation threshold for diameter polydisperse spherocylinders while 
it overestimates it for length polydisperse ones. 

Finally, we focus on connectedness distance polydispersity, illustrated in fig. [5] for the case of an 
aspect ratio of 15. As explained above, this does not lead to variability in the interactions between 
the particles, only in the definition of which particles are part of the same cluster. Nevertheless, it 
can be treated theoretically in the same way as the length and the diameter polydispersity due to the 
definition of the connectedness potetial u + . This explains why we find similar behaviour as for length 
and diameter polydispersity. Indeed, there is quasi-universal scaling with respect to A, even though 
the aspect ratio is not all that large and one would expect higher order moments in the distribution 
of connectedness ranges to show up. Indeed, our calculations show contributions up to the sixth 
moment, (£ 6 ), for the cluster size. Apparently terms involving the first and second moments, (£) and 
(£ 2 ), predominate the percolation threshold. 

Again, prediction from the second-virial approximation disagrees significantly with the simulations. 
The Lee-Parsons correction improves upon the quality of the prediction, but still overestimates the 
threshold. Both theories improve for large values of A, i.e., for globally thick connectedness shells 
for which the percolation threshold occurs at low volume fractions of particles. Our Ansatz for the 
connectedness direct correlation function has by construction the spatial structure of a second virial 
theory, even though that our Lee-Parsons extension does have the thermodynamics that goes beyond 
it. Arguably, at higher densities the actual structure of the direct correlation function starts to deviate 
from this. Note also that the Percus-Yevick prediction for the percolation threshold of monodisperse 
hard spheres is anyway in only qualitative agreement with simulations [39]. 

4.2 Cluster formation mechanisms 

As we have seen, the sensitivity of the percolation threshold to the higher order moments of the 
distribution function depends on the type of polydispersity: length, diameter or connectedness range. 
It suggests that the mechanism by which the particles cluster differs between the different types of 
polydispersity. Indeed, hard particles of different size and/or shape have for entropic reasons a tendency 
to phase separate, and even if they do not actually phase separate this might give rise to fractionation 
of particles in the transient clusters that form in the mixtures |40l ITT] . 

To investigate this, we compare in fig.[()]the distribution of lengths and diameters within the largest 
cluster V c \us{x) with that in the entire system V(x) with x = (L, D), for length and width polydisperse 
rods respectively. In both cases the larger particles are more abundant in the largest cluster than in 
the whole system, explaining why relatively small amounts of large particles have a large effect on 
the percolation threshold [§]. On the other hand, the effect weakens with increasing volume fraction 
of particles, at least for the length polydisperse ones. The proportion of short particles within large 
clusters is more and more important as packing fraction increases, making the gap between V c ius{%) 
and V(x) smaller. 

Another measure for the cluster structure is the fractal dimension df of the critical cluster. We 
obtain this by measuring the radius of gyration Rf ~ of the clusters as a function of the 
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Figure 6: Difference between the length and diameter distribution P c iuster ( x ) °f the largest cluster and the distribution 
P(x) of the whole system, normalized by the value of the distribution at its peak, P(x p ), as a function of the length 
x = L or width x = D scaled to the peak value, for two volume fractions in both cases. The length (triangles) and 

diameter (circles) distributions of the entire system are of the Weibull form P(x) oc ^ with a = 2.37 and 

/3 being such that the distribution peaks lie at L = 15 and D = 1. Larger and thicker particles cluster more easily. The 
difference between the distributions within clusters and the global ones become smaller with increasing volume fraction. 
Notice that length polydispersity has a much stronger fractionation effect than diameter polydispersity. 



Figure 7: Radius of gyration of clusters as a function of their size on a logarithmic scale for hard rods. Compared 
are results for monodisperse and for length and width polydisperse rods (see legend key). We set D = 1 and L = 15 

in the monodisperse system. Length and width distributions are of the Weibull form P(x) oc e~^ x ^^ a with 

a = 2 and (3 chosen such that {L) w = 15 and (D 2 ) = 1 in both other systems. <j> ~ (j) c and £ = 0.2 in all three 
cases. Monodisperse systems as well as systems polydisperse in diameter superimpose and can be fitted by a power law 
exhibiting an exponent around 0.35 (deviations from this scaling for very large clusters are finite-size effects). This is 
not true for length polydispersity. 
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number of particles in it, A/iuster; see fig. [3 This quantity is sensitive to length polydispersity and less 
to diameter polydispersity. The latter and that for monodisperse rods collapse exactly on the same 
curve that seems to exhibit a power law scaling and a fractal dimension of df = 2.8. Note that finite 
box-size effects cause deviations from pure power-law behaviour. For length polydisperse rods we also 
do not find power law scaling but the trend seems to conform to the same fractal dimension but with 
a larger prefactor. The fractal dimension of 2.8 is larger than the mean-field value of 2 we expect to 
hold for very long rods |421 | and that we obtain from the second virial approximation, but close to the 
accepted value of 2.5 in three dimensions for standard percolation [53]. We expect that because in 
our simulations the aspect ratio of the particles is not very large that we find a deviation from the 
mean-field exponent. 

We conclude that having varying lengths or diameters within a collection of hard rods does not 
fundamentally change the way percolation is reached. Still, the volume fraction at the percolation 
threshold is different from that for monodisperse rods, even if the average diameters (or lengths) are 
equal. 

5 Conclusions 

In summary, we have presented a theoretical and computer simulation study on the effects of poly¬ 
dispersity on the geometrical percolation in suspensions of hard spherocylinders. We compare results 
for bidisperse, Gaussian and Weilbull distributions and show that the percolation threshold is quite 
insensitive to the precise distribution. In the case of length and connectedness polydispersity the 
thresholds superpose within numerical error when scaled with the appropriate second order cumulant 
of the size distribution. For diameter polydispersity, however, higher order moments seem to matter, as 
the superposition of the different distributions is not quite perfect. To analyse the simulation results, 
we also present a theoretical treatment of the problem within connectedness percolation theory that 
we find to quantitatively predicts the percolation threshold for hard rods of aspect ratios above 20. 
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